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ABSTRACT 

We estimate the polarized thermal dust emission from MHD simulations of protostellar collapse and 
outflow formation in order to investigate alignment of outflows with magnetic fields. The polarization 
maps indicate that alignment of an outflow with the magnetic field depends on the field strength inside 
the cloud core; the direction of the outflow, projected on the plane of the sky, is aligned preferentially 
with the mean polarization vector for a cloud core with a magnetic field strength of 80 //G, while it 
does not tend to be aligned for 50 /jlG as long as the 1000 AU scale is considered. The direction of the 
magnetic field at the cloud center is probed by the direction of the outflow. In addition, the magnetic 
field at the cloud center can be revealed by ALMA even when the source is embedded deeply in 
the envelope. The Chandrasekhar-Fermi formula is examined using the polarization maps, indicating 
that the field strength predicted by the formula should be corrected by a factor of 0.24 — 0.44. The 
correction factor has a tendency to be lower for a cloud core with a weaker magnetic field. 

Subject headings: ISM: clouds — ISM: jets and outflows — MHD — stars: formation — polarization 


1. INTRODUCTION 

Magnetic fields are believed to control not only pro¬ 
tostellar collapse, but also formation of circumstellar 
disks and outflows. Many observations have suggested 
that the outflow and jet axis of the young star is 
aligned pr eferentially along the clo ud-sc ale magnetic 
Held (e.g.. ICohen. Rowland. & Blaid 11984 IStrom et all 
^22 | Vrba et Id. y~98j| IVrba 1 fftrony fo Stroml 119881: 

However, recent observations indicate a suggestion 
contrary to previous ones concerning the issue of 
alignment of outflows and jets with the magnetic 
fields. High-resolution observations of submillimeter 
polarization have resolved the magnetic fields around 
young stars on a ^ 10 3-4 AU scale, which is com- 


parable to the outflow scale dMomose et al.1 

2001 

Henniner et al.l 120011: IWolf. Launhardt. &; Hennind 

200,9 

Vallee. Greaves, fe Fie»4l2003). IWolf et alJ (2003) 

inves- 


tigate alignment of outflows with magnetic fields for Bok 
globules associated with Class 0 protostars and Class I 
sources, suggesting that two Bok globules are associated 
with outflows parallel to the magnetic fields, while the 
two other globules are associated with outflows perpen¬ 
dicular to the magnetic fields. For Cla ssical T Tauri stars 
(CTTSs), iMenard fc Duchenel ({2004 1 estimate orienta¬ 
tion of the symmetry axes of the disk-jet systems in the 
Taurus-Auriga region, and indicate that CTTSs are ori¬ 
ented rando mly with respect to the lo c al ma gnetic field. 

Recently, iMatsumotcTfe Tomisakal (120041) (hereafter 
MT04) have performed MHD simulations of the collapse 
of magnetized cloud cores and reproduced outflow gen¬ 
eration in order to investigate the directions of outflows, 
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circumstellar disk, rotation, and magnetic fields. The 
simulations show that the outflow tends to be aligned 
with a local magnetic field of a 10 AU scale irrespective 
of the magnetic field strength assumed, while the align¬ 
ment depends on the field strength on the cloud core 
scale. A disk-outflow system is aligned with the cloud 
core scale magnetic field within ~ 5° and ~ 30° for the 
initial field strengths of 37.1 and 18.6/iG, respectively, 
because of the strong magnetic braking during the col¬ 
lapse. When a weak field strength of 7.42 /iG is assumed, 
the outflow is not aligned with the cloud core scale mag¬ 
netic field. In this Letter, alignment of an outflow with 
magnetic field is discussed by constructing polarization 
maps from the MHD simulations of MT04. 

2. CLOUD MODEL AND POLARIZED EMISSION 

Polarization of dust emission is calculated from MT04 
MHD simulation data. The simulations follow the gravi¬ 
tatio nal collapse of cloud cores, formation of a first stellar 
core (lLarsonlll 969ft , and the launch of an outflow, resolv¬ 
ing both the whole cloud core and the protostar. Polar¬ 
ization maps are constructed by extracting the central 
cubic region (9128 AU) 3 from all the simulation data. 

The initial model of a cloud core is a slowly rotat¬ 
ing, spherical, isothermal cloud threaded by a uniform 
magnetic field. The initial cloud core has the density 
profil e of a Bonnor-Ebert sphere ( Ebert iBonnod 

I1956T) . and the present model can be applied to Bok glob¬ 
ules because Bok globu les are thought to have Bonnor- 
Ebert d ensity profiles dAlves et all I200H lHarvev et all 
1200 It iRacca. Gomez. & Kenvonll2002ft. The initial cen- 
tral density is set as po = 1 x 10 -19 gem -3 , which corre¬ 
sponds to a number density of no = 2.61 x 10 4 cm -3 for 
the assumed mean molecular weight of 2.3. The initial 
temperature of the gas is 10 K, the edge of the cloud is 
located at r = 0.178 pc from the center, and the mass 
of the cloud is 6.130 M 0 . The initial angular velocity is 
assumed to be 7.11 x 10 -7 yr -1 . In this Letter, the last 
stages of two models of MT04, models MF45 and WF45, 
are shown. In these models, the initial magnetic field 
is inclined at an angle of 0 = 45° from the initial rota- 
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TABLE 1 

Parameters and Properties of Model Cores 



Bo 

B 

n 

4>3D 

Models 

(mG) 

(mG) 

(cm- 3 ) 

(deg) 

MF45 

18.6 

82.8 

2.74 x 10® 

12.4 

WF45 

7.42 

50.1 

2.78 x 10® 

53.5 


tion axis, which corresponds to the z-axis. The initial 
magnetic field strengths Bo are shown in Table [l] 

Table 0 also shows the mean strength of magnetic 
field B and mean number density n within the region 
(9128 AU ) 3 for comparison with polarization maps in 
§ 03 The mean values B and n are considerably larger 
than the initial values B 0 and no because of amplification 
due to collapse. Angle (j) 3 D denotes a three-dimensional 
angle between the direction of the mean magnetic field 
within r < 50 AU and that within 4555 AU, indicating 
change in the direction of the magnetic field inside the 
cloud core (see Figs. 7 and 15 of MT04 for models MF45 
and WF45, respectively), and model WF45 shows an an¬ 
gle larger than that in model MF45 because of the weak 
magnetic braking. The direction of the mean magnetic 
held within r < 50 AU approximately coincides with the 
local direction of the outflow, and is adopted here as an 
index of the direction of the outflow. 

The polarization of du st emission is ca l culate d fol¬ 
lowing the analysi s of iFiege fc Pudritzl ((200 01 and 
IPadoan et aD (20111). As these studies do, we focus on 
the thermal dust emission at submillimeter wavelengths, 
neglecting scattering. The Stokes parameters Q and U 
are proportional to the following integrals of q and u 
when we assume that the grain properties and tempera¬ 
ture are constant: 


q = J p cos 2if cos 2 yds, 

(i) 

u= p sin 2-0 cos 2 yds, 

( 2 ) 


where the integrals f ds are performed along the line 
of sight, p denotes the gas density, if is the angle be¬ 
tween the projection of the magnetic field on the plane 
of the sky and the north, and 7 is the angle between 
the plane of the sky and the local direction of the 
magnetic field. The polarization angle x is given by 
tan2y = u/q, where — 7r/2 < x < tt/2. The polariza¬ 
tion vector with x traces the polarization shown by the 
B-vector of the submillimeter thermal emission, which 
is parallel to the interstellar magnetic field on the plane 
of the sky. The degree of polarization is calculated by 

P = a (q 2 + u 2 ) 1 / 2 / (Y, — aY, 2 ) with £ = f pds and 
£2 = (1/2) f p (cos 2 7 — 2 / 3 ) ds, where a is a parameter 
specified by the grain properties, and we adopt a co n- 
stant value of a = 0.15 following IPadoan et a 1 J ( 2001 ). 

3. RESULTS 

3.1. Polarization Maps of the Cloud Cores 

Figure [I] shows the polarization maps overlaid by the 
surface density for models MF45 and WF45 for three 
orthogonal lines of sight (x, y , and ^-directions). The el¬ 
liptical distribution of the surface density reflects the flat 


infalling envelope for all the models and for all the lines of 
sight (see Fig. 2 of MT04 for the three-dimensional den¬ 
sity and magnetic field structures of model MF45). The 
hourglass structure of magnetic fields in three dimensions 
is projected directly on the polarization vector. The ori¬ 
entations of the mean polarization vector (thick green 
line) are perpendicular to the long axes of the surface 
densities. The degree of polarization along the long axis 
is less than that along the short axis due to the hourglass 
structure of the field lines; this field structure increases 
the component of the magnetic field parallel to the line of 
sight Bios compared with the perpendicular component 
B ±, and reduces the degree of polarization along the flat 
envelope. 

The outflows are hardly visible in the maps because 
the outflows are extended up to only ^ 200 and 150 AU 
for models MF45 and WF45. We adopt here the direc¬ 
tion of the mean magnetic field within r < 50 AU as 
an index of the direction of the outflow (thick red line). 
The projected direction of the outflow is almost aligned 
with the mean polarization vector (thick green line) for 
model MF45 irrespective of the lines of sight (Fig. [Ha- c) 
because of a small intrinsic angle between the outflow 
and the cloud core scale magnetic field (03 d = 12.4°). 
On the other hand, the alignment depends on the line- 
of-sight for model WF45 (Fig. EJT/) because the out¬ 
flow is not aligned with the magnetic field intrincically 
(03D = 53.5°). When model WF45 is observed along 
the x-direction (Fig. [H^), the projected direction of the 
outflow is aligned with the mean polarization vector by 
chance. They are not aligned considerably when ob¬ 
served along the y- and ^-directions as shown in Figure's 
and^Jf 

More qualitatively, Table0shows 0 pr j, which is defined 
as an angle between the red and thick green lines. This 
angle indicates that the direction of the outflow tends 
to be aligned with the mean polarization vector on the 
plane of the sky when a cloud core has a strong magnetic 
field of B ~ 80/iG. 

The polarization maps examined here demonstrate 
that the direction of the magnetic field at the cloud center 
can not be inferred from the polarization on the scale of 
this map (^ 10 3 AU scale). Model MF45 shows polariza¬ 
tion maps similar to model WF45, although the magnetic 
field at the cloud center of Model MF45 is projected to 
quite different direction from that of Model WF45.The 
directions of the central magnetic field reflect the those 
of outflows. In other words, the direction of the outflow 
probes the direction of the magnetic field at the cloud 
center. 

Figure [ 2 simulates direct observations toward the cloud 
center on the 100 AU scale with high resolution of 25 AU 
for the three lines of sight for model WF45. These 
maps are constructed by convolution of the polariza¬ 
tion pattern of Figure EETEEK with a Gaussian beam of 
FWHM = 25 AU, which corresponds to angular resolu¬ 
tion of 0.1" for an object at a distance of 250 pc, e.g., an 
observation toward B335 by ALMA. In constructing Fig- 
ureO the envelope of 9128 AU scale is taken into account 
similar to Figure On the scale of Figure 0 the po¬ 

larization reveals the convergence of alignment between 
the disk-outflow system and the magnetic field; never¬ 
theless the protostar is embedded deeply in the envelope. 
Figure 0 shows the polarization depending on the lines of 
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Fig. 1. — Polarization maps constructed from MHD data of the central (9182 AU) 3 box for models (a-c) MF45 and ( d-f) WF45 along 
the lines of sight parallel to (a, d) the cc-axis, (b, e) the y- axis, (c, f) and the z-axis, respectively. Vector traces the B-vector of the polarized 
thermal emission, and is therefore parallel to the interstellar magnetic field lines on the plane of the sky. The length of the vector is 
proportional to the degree of polarization, which is also shown by black contours at 1 % intervals. The mean orientation of the polarization 
vector is denoted by the thick green line. Color scale and blue contours denote the surface density E, which is proportional to the intensity 
of thermal dust emission when the gas is assumed to be optically thin. Thick red line denotes the projected direction of the mean magnetic 
fields averaged over r < 50AU, indicating the direction of the outflow. 


TABLE 2 

Outflow Orientations and Estimate of Magnetic Field Strengths 


Models 

Line of 

Sight 

A'.i 

(deg) 

P 

(g cm -3 ) 

Sv 

(km s _1 ) 

O'x 

(deg) 

Bcf 

( m g) 

/ 

MF45 

X 

4.40 

1.04 x 10- 18 

0.209 

10.3 

242.8 

0.34 

MF45 

y 

13.1 

1.04 x 10- 18 

0.214 

11.0 

233.2 

0.35 

MF45 

z 

9.62 

1.04 x 10- 18 

0.193 

12.4 

186.4 

0.44 

WF45 

X 

-3.02 

1.06 x 10- 18 

0.262 

16.0 

197.6 

0.25 

WF45 

y 

59.1 

1.06 x 10- 18 

0.251 

14.3 

212.2 

0.24 

WF45 

z 

86.5 

1.06 x lO" 18 

0.244 

17.2 

171.6 

0.29 


sight (see also Fig. 13 of MT04 for the three-dimensional 
structure). Figured shows an edge-on view of the disk 
as reproducing the considerably flat surface density, and 
the polarization vector is almost perpendicular to the 
long axis of the surface density. The bipolar outflow is 
extended up to 150 AU vertically on the map, although it 
is hardly visible in both the polarization and the surface 
density. Figured also shows an edge-on view, exhibit¬ 
ing a polarization pattern different from Figure EJa; the 
polarization vector is oriented along the projected disk 
surface because the radial component of the magnetic 
field in the hourglass structure contributes B± there. In 
\z\ > 100 AU, the polarization vector is still perpendicu¬ 


lar to the disk, similar to Figure Eb- Figure (2t shows the 
change in the direction of the mean magnetic field; the 
outflow, whose direction is indicated by the thick line, is 
aligned with the central polarization vector. 


3.2. Estimate of Field Strengths from Polarization Maps 

Magnetic field strength can be estimated from the 
polarization maps of F i gure ED using the method of 
iChandrasekhar fc Fermi 


£?cf 



(3) 


log E [cm" 
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Fig. 2.— Polarization maps convoluted by the Gaussian beam with FWHM = 25AU toward the cloud center for model WF45 along 
the lines of sight parallel to (a) the cc-axis, ( b ) the y-axis, and (c) the z-axis. The beam pattern is denoted in the bottom left corner of 
each panel. Thick line denotes the projected direction of the mean magnetic fields averaged over r < 50AU, indicating the direction of the 
outflow. Grayscale and contours denote the surface density. 


where p denotes the mean density, estimated by averag¬ 
ing the surface density over the map, Sv denotes the rms 
gas velocity, assumed as the velocity component parallel 
to the line of sight superposed as Sv 2 = f pv 2 os dV / f pdV, 
and a x denotes the standard deviation to the mean ori¬ 
entation angle of the polarization vector. The derived 
parameters p, Sv, a x , and the estimated magnetic field 
strength Bqf are shown in Table O for the three orthog¬ 
onal lines of sight. 

The estimated field strengths range from 186 to 242 pG 
for model MF45, and from 171 to 212 pG for model 
WF45, exhibiting significant dispersion 5 . Moreover, 
these field strengths are 2 — 4 times larger than the mean 
magnetic field strengths B (see Tabled, which are ob¬ 
tained directly from MHD data. In other words, the 
prediction of the Chandrasekhar-Fermi formula should 
be corrected by a factor of / = 0.34 — 0.44 for model 
MF45 and / — 0.24 — 0. 29 for model WF45, as s hown 
in Table □ Jpstffiey^pne, & Cammii and 

iPadoan et al.1 ( 200 11 also report such an overestimated 
field strength in application of the Chandrasekhar-Fermi 
formula, and obtain a correction factor of f ~ 0.4 — 
0.5, which is consistent with our case. IPadoan et al.l 
( 20011 discuss the deviation from the prediction of the 
Chandrasekhar-Fermi formula as being attributable to 
deviation from the linear theo ry, which is assumed by 
iChandrasekhar fc Fermil (11953 1. The large values of a x 
shown in Table El indicate that the Alfven wave is non¬ 
linear with large amplitude in our models. Moreover, it 
is noteworthy that model WF45 exhibits larger a x than 
model MF45, indicating that the correction is more sig¬ 
nificant for a cloud with weaker magnetic field. 

4. DISCUSSION 

Activity of an outflow may be discussed in terms of 
alignment of the outflow with the magnetic field. The 
cloud core with stronger magnetic field exhibits a faster 
outflow as shown in Figure 18 of MT04, and the out¬ 
flow tends to be aligned with the polarization vector. 

5 We confirmed that a method proposed by IHoudd <2001 can 
correct a projection effect of the magnetic field. 


This indicates that the fast outflow is observed parallel 
to the magneticyfield. According to the o bservations to¬ 
ward CTTSs, iMenard fc Duchenel (120041 1 suggest a sim¬ 
ilar tendency in spite of the different evolutionary stage 
from that considered here: CTTSs without bright and 
extended outflows have a tendency to be perpendicular 
to the magnetic field. 

The outflows are extended up to only 150 — 200 AU in 
the MHD simulation data used here, while the outflows 
observed by molecular line emission are extended up to 
a 1000 AU scale (e.g.. lWolf et alJl2003 f. Moreover, the 
outflows presented here have considerably slower veloc¬ 
ity than that observed around young stars. Therefore the 
cloud cores presented here are restricted to the very early 
evolutionary stage compared with the observed cloud 
cores. This restriction arises in response to computa¬ 
tional cost for solving the launch mechanism of the out¬ 
flow near the protostar. In further stages, the magnetic 
field hardly seems to affect the direction of the outflow 
on the scale of 1000 AU, because the outflow is acceler¬ 
ated at r ~ 10 AU at a speed comparable to the Alfven 
velocity at this radius, and the magnetic field strength of 
the envelope decreases rapidly, proportional toFocr -1 . 
Moreover, the protostar may be decoupled from the mag¬ 
netic field of the envelope as a result of efficient ambipolar 
diffusion iXakano. Nishi. fc UmebavashifiOQ^ in further 
stages. The directions of the outflow will be fixed during 
the main accretion phase of the protostar. 

Numerical computations were carried out on the 
VPP5000 supercomputer at the Astronomical Data 
Analysis Center (ADAC) of the National Astronomical 
Observatory, Japan. This research was supported in 
part by Grants-in-Aid for Young Scientists (B) 16740115 
(TM), for Scientific Research (C) 14540233 (KT) and 
17540212 (TM), and for Scientific Research (B) 17340059 
(TM, KT) by the Ministry of Education, Culture, Sports, 
Science and Technology, Japan. 
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